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In this paper, we give a very accurate description of the way the 
simple exclusion process relaxes to equilibrium. Let Pt denote the 
semi-group associated the exclusion on the circle with 2N sites and 
N particles. For any initial condition x, and for any t > log N, we 
show that the probability density Pt(x, •) is given by an exponential 
tilt of the equilibrium measure by the main eigenfunction of the par¬ 
ticle system. As log N is smaller than the mixing time which is 
logN, this allows to give a sharp description of the cutoff profile: 
if djv(t) denote the total-variation distance starting from the worse 
initial condition we have 

lim djvf^logiV-h^s') =erf('22g-A 
N^oo J \ n J 

where erf is the Gauss error function. 

1. Introduction. The exclusion process is a lattice interacting particle 
system where particles perform independent nearest-neighbor random walks 
with the added constraint that each site can be occupied by at most one 
particle (see the classic references [13] and [14] for a complete introduction to 
the subject). It is a very simplified model for a gas of interacting particles. 
We consider in this Introduction the case were the lattice is either a d- 
dimensional torus or hypercube of side length N. The number of particle is 
chosen to be proportional to the number of sites. 

In this paper, we investigate how the particle system starting far away 
from equilibrium, relaxes to its equilibrium state. This question can in fact 
be treated with different point of views: 

• One can describe the evolution of the system at the macroscopic level, 
giving the evolution of the density of particle in space and time. This 
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is the study of hydrodynamic limits (see [7] for an introduction to the 
subject). 

• One can adopt a microscopic point of view, and look at the evolution of 
the law of the particle, and in particular, its total variation distance to 
the equilibrium law. This is the study of the Markov chain’s mixing time 
(see [12]). 

With the first point of view, it is now well known that under diffusive rescal¬ 
ing (space rescaled by N and time rescaled buy N'^), the density of particle 
evolves like the solution of the heat equation. The result is valid in any 
dimension (see [7] for references). 

Concerning the mixing-time approach, progresses are more recent. It has 
been shown by Morris that in any dimension the time needed to come close 
to equilibrium in total variation was of order A^^logA^ [16]. In dimension 1, 
more refined estimates have been obtained and gave the exact location of 
the mixing time either for the segment [8] or the circle [9] with lower bounds 
proved earlier by Wilson [19] (see also [17] for results in the case of arbitrary 
graph, and [6, 10]). A natural question is then what does the law of the 
particle system look like when it is about to reach equilibrium. 

At equilibrium, the law of the distribution is uniform over all particle 
configurations. Another way to see it is to say that the state of each site 
(occupied or vacant) is given by a field of i.i.d. Bernoulli variables whose 
sum is conditioned to be equal to the number of particle (which is a fixed 
parameter). 

What would be natural to expect then, is that up to a small correction, 
before equilibrium, the particle distribution still is a conditioned product 
measure, but that the Bernoulli variables are no more identically distributed: 
there is a space dependent bias which is given by the solution of the heat 
equation. This brings a strong connection between the problem of the mixing 
time and that of the hydrodynamical limit. This connection was previously 
underlined by Lee and Yau when studying the related issue of log-Sobolev 
constant for the simple exclusion [11]. Indeed in the case of small bias, with 
some minor efforts one can derive sharp estimates on the total-variation 
distance between the conditioned product of biased Bernoulli and the equi¬ 
librium measure. This can be turned into a precise prediction on how the 
total-variation distance drops from one to zero, the cutoff profile. The present 
paper brings this heuristic picture on a rigorous ground in the case of the 
exclusion on the circle. 

2. Model and results. 

2.1. The process. We consider Z 2 Ar := Z/(2iVZ), the discrete circle with 
2N sites and we place N particles on it, with at most one particle per 
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site. With a slight abuse of notation, we will sometimes use elements of 
2N} C Z to refer to elements of Z 2 n- 
The exclusion process on 7j2n is a dynamical evolution of the particle 
system which can be described informally as follows: each particle tries 
to jump independently on its neighbors with transition rates p{x,x + 1 ) = 
p{x,x — 1 ) = 1 , but the jumps are canceled if a particle tries to jump on a 
site which is already occupied. 

Let us describe the chain more formally. We adopt the convention that 1 
denotes a particle and —1 denotes an empty site. This is not the most usual 
one (empty sites are more often denoted by 0 ) but it proves to be more 
practical in our computations. Our state-space is dehned by 

{ 2N 

X=1 

Given 77 G 0 define the conhguration obtained by exchanging the content 
of site X and x + 1 

{ rj^{x) := r]{x + 1), 
rj^lx + 1) := r]{x), 
r]^(y)=il(y) Vy^{x,x + 1}. 

The exclusion process on 'Z 2 N with N particle is the continuous time Markov 
process on Otv whose generator is given by 

(2-3) {CNf){r]):= ^ /(??"')-/(??)• 

xeZ2N 

The chain is irreducible and reversible, and the unique invariant probability 
measure is the uniform probability measure on ftjsf which we denote hy pisf. 
Given x we let {r]f)t>o denote the trajectory of the Markov chain 

starting from x- We write P[(r 7 ^)t>o £ •] for the law of (? 7 ^)t>o- We let Pt 
denote the Markov semi-group and write for the probability measure 

Ptixr), X^^N- 

We measure the distance to equilibrium in terms of total variation dis¬ 
tance. If a and /3 are two probability measures on fl, the total variation 
distance between a and f3 is dehned to be 

(2.4) ||a - /3 ||tv := ^ ~ 

where x+ = max(x, 0) is the positive part of x. It measures how well one can 
couple two variables with law a and 13. We dehne the distance to equilibrium 
of the Markov chain to be 


d^it) := 


maxIlP/^ -//|| 

TV 


(2.5) 
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In a previous paper [9], we described in detail the asymptotic behavior 
of d^{t). We proved that around a time of order ^logA^ the distance to 
equilibrium drops from 1 to 0 in a time window of width N'^, 

lim limsupd^ f log =0, 

5^.00 \^27r^ J 

2.6 

/ Ar2 \ 

lim liminfd^l —^logA^ —A^^s) =1. 

s—>—oo N^oo V 27r2 / 


The aim of this paper is to complete this picture by identifying, in an 
acute way, the pattern of relaxation to equilibrium. In particular, we are 
interested in proving the existence and finding an expression for limiting 
profile 

(2.7) lim { —^-logA^ + A^^s 

N^oo \2 tX^ 

To reach this goal, we have to understand what the distribution looks 
like much before the time ^^logA^. 


2.2. The mixing time profile. The main achievement of our paper is to 
determine the cutoff profile. 


Theorem 2.1. The total-variation distance to equilibrium from the worst 
initial condition has the following asymptotic profile: for any s G M, we have 

(2.8) lim d^(^\ogN + ^s\= erf I , 

where erf is the Gauss error function 

(2.9) erf(t) dn. 

A Jo 

The method by which we obtain the result gives us in fact much more 
information about the relaxation of the system: we are able to characterize 
fully how Pfi looks like much before the mixing time, for all initial condition 
X G Qn- 


Remark 2.2. The fact that the prohle of the cutoff is given by a function 
of the type erf(^e“*) (the constant is not essential since it depends on 
the particular process and the choice for the normalization) is given by 
Wilson [19] as a conjecture (supported by numerical evidences) for a process 
very much related to the exclusion: the adjacent transposition shuffle. The 
reason why the function erf appears is that the last statistic that comes 
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to equilibrium for the process (here the first Fourier coefficient of rj, see 
below) is well approximated by a Gaussian; the exponential terms are present 
because the mean of this Gaussian converges exponentially to zero. This 
is a property which is believed to be shared by many Markov chains and 
rigorously known, for example, the random walk on the hypercube [5]. Let us 
mention that however there are the Markov chains which exhibit cutoff and 
do not have this property. This is, for instance, the case of top to random 
shuffle [4], and also of the transposition shuffle for which the important 
statistic, the number of fixed point, behaves like a Poisson variable (see, 
e.g., [15]). 

2.3. The description of Pf" much before equilibrium. The main result of 
the paper, from which we deduce Theorem 2.1 requires some notation to 
be introduced. The time evolution of the density of particles is given by 
the discrete heat equation on 'Z 2 N and for this reason, the eigenfunction of 
the discrete Laplacian on the circle plays an important role in our analysis; 
in particular, those in the eigenspace corresponding to the spectral gap: 
X i-T- cos(^), and x sin(^). 

To describe the projection of y G 11^ on this eigenspace, it is more conve¬ 
nient for us to have one positive coefficient than two real ones, and for this 
reason we introduce 0{x) which is the “phase” of y in the first eigenspace. 
It is the unique 0 G [0,27r) satisfying 



or 6{x) = 0 if the system has no solution. We denote by b{x) the hrst Fourier 
coefficient of x 

(2-11) 6(x):=^ x{x)sm(^^ + 9 

zG'Ia2N 

Note that b(x) = 0 in the case where (2.10) has no solution. In the case 
where x = +1 for x G {1,... ,N} and —1 elsewhere, r]{x) = 

If ^ is a probability measure on a state-space 11 and that f is a function 
H —)• M, we use the following notation for the expectation: 

v{f) ■■= ^{fih)) ■■= fivXh)- 


( 2 . 12 ) 


6 


H. LACOIN 


We define given a > 0 and 6 G [0,27r). We define to be the proba¬ 

bility measure on with the following Radon-Nikodym density: 


(2.13) 
where 

(2.14) 


Finally, let us set 
(2.15) 


^j^N,a,e ^aaeirj) 

duN /iAr(e"“9(^))’ 


ae(??) := y~]7?(x)sin 



A^:=2(l-cos(^)). 


We note that Atv is the spectral gap of the simple random walk on 'Z 2 N 
(with jump rate one in each direction), and hence from [3], Section 4.1.1, it 
is also the spectral gap of the exclusion process on 1 i 2 n- 

The main result of the paper tells us that much before the mixing time, 
is close to for an appropriate choice of a and 9. 


Proposition 2.3. For all N sufficiently large, for all x £ for all 
t > '^6 have 

(2.16) \\Pf - < (loglogA^)-^ 

Theorem 2.1 follows from Proposition 2.3 by using the following lemma. 


Lemma 2.4. For all K > 0, for all N sufficiently large (depending on 

K), 


(2.17) 


lim sup 
N^ao ^^[o^K] 

0e[O,27r) 




IiaIItv 



= 0 . 


2.4. Exclusion with an arbitrary number of particle. We have chosen to 
present here the result only in the case where the number of particles is 
equal to half of the number of sites. However, mutatis mutandis, the proof 
adapts to the case of particle k]\f < N on Z 27 V where /ctv tends to infinity 
with N (the case k> N can be treated by symmetry). Let us discuss here 
what the results are in that case and how they can be obtained. 

While the case of proportional to N can be derived directly from 
the proof presented here, it turns out that some of the technical lemmas 
(e.g.. Proposition 6.3) breaks down if kjy grows much slower, that is, like 
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logA^. However, adapting the techniques developed specifically for the case 
of slowly growing in [9], all technical obstacles can be overcome. 

To close this discussion, let us mention what the cutoff profiles are in 
those cases. When = [aiV] for some a = (0,1), we have 

sin(a7r/2) 


(2.18) lim d^{-^logN-\ - 5 -s)=erf 

N-^oo yZTT^ TT^ J 

When k^ satisfies 1 <C /cjv ^ we have 

iV 2 ^2 

d^' I 

A^—)-oo 




(2.19) 


lim d^ ( 7 —^ log kjsf H- 5 -s ) = erf ( -e 




vr^ 


2.5. Organization of the paper. In Section 3, we prove Theorem 2.1 from 
Proposition 2.3, and also give a proof of Lemma 2.4. In Section 4, we decom¬ 
pose the proof of Proposition 2.3 into three key statements, whose proofs 
are, respectively, given in Sections 5, 6 and 7. 


3. The cutoff profile. 

3.1. Proof of Theorem 2.1. 
check that for 


Let s G . 
Ar 2 


be fixed. It is straightforward to 


Ar 2 

ls,N := logiVH- 5 -s 


we have 
(3.1) 


sup sup 
N xeUjv 


y/Nb{x)e < 00. 


Hence, using the triangular inequalities. Proposition 2.3 and Lemma 2.4 we 
have for all x S ^N, 

71 


(3.2) 


lim 




TV 


— erf 




= 0 . 


The asymptotic for d^^tg^i^) follows if one can identify x which maximizes 
b{x)- A few seconds of thoughts are enough to realize that Xmax defined as 
follows is the unique maximizer up to translation: 

-|- 1 , for x = 1 ,..., A^, 

— 1 , for x = A"-L 1 ,..., 2 A^. 

The asymptotic behavior of l)(xmax) is given by the following limit; 

' XTT 


(3.3) 


Xmax(T) = 


(3.4) 


N E Xmax( 


N—>CXD N 

x(PIj2n 
we have also 


X sm 


N 


TT \ 4 


TT 


As for any s G ’ 

(3.5) 

the result follows from the continuity of the error function. 


lim y/Ne ^ ^ s 

N—^oo 
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3.2. Proof of Lemma 2.f. The underlying idea is quite simple: we want 
to prove that asymptotically under /tat, once rescaled 

converges to a Gaussian. 


Lemma 3.1. The following statements hold true: 

(i) For a fixed 6 G [0,27r). The quantity converges in law to a 

standard Gaussian. Moreover, the convergence is uniform in 6, in the sense 
that for any bounded continuous function F the convergence 


(3.6) 


lim uat 

N^oo 


F 


ae{v) \ 

VnJ 




/^du, 


holds uniformly in 9. 

(ii) Moreover, ag{r]) is exponentially concentrated in the sense that there 
exists a constant c > 0 such that 


(3.7) 


LN{\ad{r])\>VNu) <2e 


Let us explain how we deduce Lemma 2.4 from these facts. We note that 


(3.8) 


A,7A-i/2,0 




- 1 


Because of convergence in probability and exponential tightness, we have 

(3.9) lim = 

N^oo 

Thus, — TnWtv converges uniformly in 7 G [0,K] and in 0, to 

(3.10) 

The conclusion then follows by performing a few changes of variables. 


1 f |e7«-7V2 _ 


Proof of Lemma 3.1. Statement (ii) follows from a more general state¬ 
ment on concentration for Lipshitz functional on Lljsf. Proposition B.l is 
proved in the Appendix. 

To ensure that the convergence holds uniformly in 9, the reader can check 
that all the bounds present in the proof do not depend on 9. In the remainder 
of the paper, we will use the notation 
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As most computations do not depend on the value of 6, we omit it from the 
notation most of the time in the remainder of the paper. Note that a{r]) sat¬ 
isfies trivially fiNiaiv)) = 0 - Let us show that the variance is asymptotically 
equivalent to N. 

(3.12) = E sin(x)^-|- sin(x)sin(y)E[ry(a:)r 7 (y)]. 

X&Z2N {x,y)GZ2N 

x^y 

The first term is equal to N. As for the second term, as we have ¥^[r](x)r](y)] = 
— 1/(2A^ + 1) it is equal to 

(3.13) —-- 7 > sin(x)^ = —--. 

^ ^ ( 2 Ar- 1 ) ^ ^ 2 Ai -1 

' xeZ2N 

To show the convergence to a Gaussian variable, we will use the martin¬ 
gale central limit theorem [2], Theorem 1. Let (Afj^)^^~^ be the martingale 
dehned by 

(3.14) :=hN(a(h)l(h(a;))L=i)- 

It satishes = 0 and = a{rj). Set 

(3.15) AMi := 
and 

(3.16) al^ = yN{{AMif\{r]{x))l^^). 

To apply the central limit theorem the martingale Mi, one must only check 
that 

2N-2 

(3-17) Jv ^ 

i=0 

converges to one, in probability (there are in fact other assumptions to check; 
see [ 2 ] but in our case they are trivially satisfied). 

For A C we let r]{A) denote the number of particles in the set A, 

(3.18) r?(A):= 

xGA 


Let us fix i and set k = k{y,i) := ? 7 ([l,z]). A simple computation gives 

2N 

[X] 


(3.19) AMi = < 


N + k — i / —,. ,. 

— sm(z -F 1) — 


2N-i 

if r]{i -F 1) = 1, 




2N-i- 


- - sin(a 

-i-l ^ 


N-k 1 

2jvV7 + 


x=i+2 

2N 


X] sin(x) , 


if r]{i -F 1) = —1. 


x=i-\-2 
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As the first and second option in (3.18) have respective probability {N — 
k) / (2N — i) and {N + k — i)/{2N — i), we have 

2 + k-i){N-k) 


(3.20) = 


{2N -if 


sin(z + 1) — 


2N-i- 


2N 

7—J E 


x=i+2 

Once this is done, we just need to check the following facts to conclude: 

(a) For all i, is almost surely smaller than 8. 

(b) For all N sufficiently large, for all i G [0, 2N — VN] we have 

2 \ 




(3.21) 


1 


^LN - o sin(i + 1) - 


< 1V-V20_ 


1 


2N 


2N-i-l 


sin(x) 


x=i-\-2 


> Ar-V20 


(c) We have the following convergence: 

2N-1 / 2N 

(3.22) lim — y iffi(z + 1) - —-:-- V 

A —>00 N \ 2A^ — i — 1 

i=l \ x=i-\-2 

From the these three claims, it is rather standard to show that converges 
to 1 in probability and we leave it as an exercise to the reader. Item (a) is 
obvious, item (b) follows from computing the mean and variance of k{r],i) 
[which are, resp., equal to i/2 and i(2A^ — i)/4(2A^ — 1)] and applying the 
Markov inequality. As for (c), it can be checked via a tedious but simple 
computation. □ 



4. Decomposing the proof of Proposition 2.3. We present in this section 
the main steps of the proof of our main result and the heuristics behind it. 


4.1. Why coupling with ? First, let us try to understand why 1 /^’°^’^ 

gives a good approximation of the Pf , via an analysis of the particle density 
and fluctuation. Let 


(4.1) u^{x,t) :=E[t]/{x)] 

denote the expected particle density (with our convention it can be negative 
since empty sites count for —1). It is standard to check that is the solution 
of the discrete heat-equation 

. . (dtu{x,t) := Au{x,t), 

^ ’ \uix,0):=xix), 

where A denotes the discrete Laplacian 

(4.3) Au{x,t) := u{x -|- l,t) + u{x — l,t) — 2u{x,t). 
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Here and in what follows if / is a function of 'L 2 N (identified to a periodic 
function of Z) such that 

(4.4) ^ /(t) = 0 , 

X^'L2N 

and X and y are two elements of 'L 2 N and xq < yo two elements of Z which 
are, respectively, equal to x and y modulo 2N. Then we use the notation 
denote the sum fi^)- straightforward that it does 

not depend on the particular choice of xq and yo once x and y are hxed. Let 
us quickly investigate the fluctuations of the integrated density of particle 

X 

(4.5) ^{t]){x) :=J2v{z)- 

Z = 1 

At equilibrium, ^{r]) is a simple random-walk conditioned to return to zero 
after 2N steps. It is centered and has Gaussian fluctuations of order ^/N. In 
[9], we have proved that the fluctuation of i{r]t){x) around its mean [given 
by ™ always of order y/N. 

This gives the intuition that much before the mixing time, the law of 
r]f might approximately be that of 2N independent ±1 Bernoulli variables, 
each with bias u^{x,t), conditioned to 

For t > log A^, u^{x) is very well approximated by a sinusoid function 
(see Lemma 5.3) 

(4.6) u^{x, t) « b{x)e~^^* ’ 

and the conditioned law of independent Bernoulli described above is very 
close in total variation to _ 


4.2. How to do it. Let us first write here the rigorous result concerning 
the fluctuation around the expected density of particle. 

Proposition 4.1. There exists a constant c > 0 such that for all N 
sufficiently large, for all t > |(AAr)~^ log AI, we have 


P? 


(4.7) 


3x,y G Z2Ar, 


<2e“ 


E 


z=x-\-l 


r]{z)-e ^^^b{x)sm{ — + 9{x) 


> sVn 
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In particular, we know that with high probability, rjf lies in the set 




max 


(4.8) 


-asin( ^ + 0 


Z=X-\-l ^ 

< Vn log log N I 






with a and 0 being chosen, respectively, equal to e 9{x)- 

To prove Proposition 2.3, it is sufficient to prove that: 

• within a time A^^(log (i-e., a shorter time-scale than the mixing 

time), one can couple a dynamics with initial condition y S with one 
with initial condition distributed like 

• the family of measure is almost preserved by the dynamics in the 

sense that applying the semi-group Pt to it only changes the value of a. 


Both of these statements hold provided a is sufficiently small, and are stated 
as two propositions below. More precisely Let be the law of a system 

started with initial conhguration 

(4.9) z/f’"’'^(ry) := ^ r-f’“’^Pz(7 /',t/). 


Proposition 4.2. For all N sufficiently large, for all 6 G [0,27r), for all 
a < , we have for all x G Ga,9, for all t > N'^{logN)^/‘^ 

(4.10) ^ 


2 log log N 


Proposition 4.3. There exists a constant C such that for all N and 
all a>0, 

//I 1 1 \ II N,a,9 N,ae~^n^ 9\\ ^ 2 fXr 

(4.11) SUp||l/j ’ ’ —U’ ’ ||rpy < Ca V A. 

t>o 


Proof of Proposition 2.3. We have for any x in /S for t > := 

3/7(AAr)~^ we have, using the triangular inequality 


(4.12) 


\\P?- 


NMx)e-^N\9(x)\ 


I TV 


x'efijv 

-.N 


TV 


< ^ ^6(x)Af-3/r,0(x)) 

+ , max ||P£. - 


,NMx)e-^N\e{x)\ 


x'ee^ ,,, 
6(x)iV-3/7,e(x) 


ITV 
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According to Proposition 4.1, we have 

We note that for x' G ^b{x)N-^/’^,e{x) 


ipx' _,.N,b{x)e-^N\e{x)\\ 

IW-io 11 TV 


(4.14) 


< II PX _ y 


iV,b(x)N-3/h0(x)| 


I TV 


+ 


Af,b(x)N-3/’',6»(x) _ ,.N,b(x)e-^N\e{x)\ 


ITV 


The hrst term is smaller than 2 iogiogN according to Proposition 4.2 as soon 


as 


t > to + x/^ogN. 

Proposition 4.3 ensures that the second term is smaller than (log hence 

we can conclude. □ 


5. Proof of Proposition 4.1. This statement is in fact mostly derived 
from the statement about fluctuations proved in [9] which we state now. 

Proposition 5.1 ([9], Proposition 4.1). There exists a constant oO 
such that for all t>0, for all s>0, for all x G we have 


(5.1) 

3x,y£ I, 2 n, 

i'n{z)-vF{z,t)) 

IV 

v 



z=x-\-l 



Remark 5.2. Note that in [9], t > 3N‘^ is required (that would be in 
fact t > 12N'^ in our setup because we work on Z 2 Ar instead of Zat), but 
this is only to treat the case of an arbitrary number of particles. The reader 
can check from the proof that this assumption is only needed to check [9], 
equation (4.4), which is obviously valid for all t > 0 when we have N particles 
on 2N sites. 


With this result. Proposition 4.1 follows from a basic analysis of the 
Fourier coefficients of the solution of (4.2). 


Lemma 5.3. For all t > |(AAr) ^ log A^, we have 


(5.2) 


max 

X&,2N 


u^{x,t)-e +e{x) 


<4iV-V2. 
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Proof. Let us decompose along its Fourier modes for the heat- 
equation. 

As in Section 2.3, we prefer to have only one coefficient per eigenspace, and 
thus, for i = 2 ,..., — 1, introduce Oi(x) to be the phase of the projection. 

We let Oi{x) be either the unique solution of 


x(x)cos 


(i 


- 


ITTX 


(5.3) 


x£li2N 


V N 


+ ^ — 0 , 


^ x(x)sin( ^ + 0) >0 


xG'^2N 

or 0 if all 6 are solution, and let us set 


1 


X{x)sm\^ — +9i{x) 


(5.4) 

We have 

(5.5) 


z£'E2N 


ITTX 


N 




Z&2N 


N-1 


i=l 


xix) ■■= ii(x)sin( — + 0iix) ) +^iv(x)(-l)'"''- 


As the functions x i-A sin(^ -|- 9i) are eigenfunctions of the discrete Lapla- 
cian with respective eigenvalues. 


(5.6) 

we have for all t > 0, 

(5.7) 
where 


- Xi^N := 2( 1 - cos I — 


ZTT 


N 


TTX 


u{x,t) = e ^ 6(x)sin( — + 9{x) ] +R{x,t,x), 


N-l 


(5.8) R{x,t,x) ■=Y^ ^'’’^%ix) 


i=2 


f ilTX 

'\l^ 


+ ^*(x) ) + ^7v(x)e 


Noticing that |5j(x)| < 2 and that for all N >2, 

(5.9) Vi G {2,..., AT}, Xi,N>iXN, 

we have for all t > (Aat)"^ 


N 


(5.10) 


\R{x,t,x)\ < 2^( 


,— iX]\^t _ 


2e 


— 2X]\^t 


i=2 


\—Q—^Ni 


<4e 


—2A f<[t 


Hence, we have the result. □ 
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6. Proof of Proposition 4.2. We assume without loss of generality that 
a is nonnegative, and write v°^ and vf for and and a(r/) for 

aoiv)- 

6.1. Properties of . In this section, we check several properties for 
While the results are quite intuitive, their proof is quite technical and we 
have decided to postpone them to Appendix A. First, we want to ensure 
that it has the right density of the particle. 

Proposition 6.1. There exists a constant C such that for all a <1 we 
have 

(6.1) sup \u°‘{r](x)) — asm{x)\ < 0 ( 0 ^ + N~‘^). 

xeZ2N 


Then we have to check that the fluctuations are not larger than ^/N. 

Proposition 6.2. There exists constant c such that for all N > 0, for 
all \a\ < , and t > 0 


(6.2) i/f 

3x,ye I, 2 n, 

y 

{rj{z) — ae~^’^^sm{z)) 

> sVn 



z=x-\-l 



Finally, we want to check that if one starts from distribution there is 
a positive density of sites where 77(2:) / Tf(z + 1), that is, of locations where 
jumps of the particle can occur. The utility of such a statement will be 
become clear in the next section when we construct the dynamical coupling. 
For a probability measure p defined on fljy, we let be the law of the 
Markov chain (r/t)t>o starting from r/o distributed like n. Set 

(6.3) j{x,y,v) := {z e [x,y]\T]{z) r]{z + 1)} 
and 

(6.4) £ := {rjenN\y{x,y) £ Z 2 N,#[x,y] > ^ j{x,y,ri) > lif[x,y]}. 

Proposition 6.3. There exist a constant c such that for N sufficiently 
large, for all |a| < cN~^/^ 

(6.5) 

Remark 6.4. The power exponents for N in Proposition 6.3 are rather 
arbitrary and other choices would also fit. The important result is that the 
probability tends to zero. 
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6.2. The ^ dynamics. We introduce in this section an auxiliary dynamics 
(the same as in [9]) which is used to couple with x^Oa (we use this 
notation for Q^q) with uf. The idea of using interface dynamics to study 
particle system dates is not new and is already present in the seminal paper 
of Rost about the asymmetric exclusion on the line [18] (for the use of 
this technique for mixing time related issues, see [8, 9, 19]). In [8, 19], the 
height function is introduced mainly to have a better intuition on an order 
which can be defined without the interface representation. Let us stress 
that here, on the contrary, the interface dynamics is used to perform a 
monotone coupling that could not be constructed by considering only the 
original chain. 

Let us consider the set of discrete height functions of the circle. 

(6.6) rijy := ^ ^|?(0) G 2Z, Vx £ 'Tj2Ni |'^(3^) ~ T 1)1 ~ 1}- 


Given ^ in we dehne as 


f ?"'(y) = C(y), 

1 = i{x + 1) +i{x - 1) - 2^(x). 


We let be the irreducible Markov chain on whose transition rates p 
are given by 


{p{U^) = i, yxe'LN, 

\p(c,n = o, 


We call this dynamics the corner-flip dynamics, as the transition ^ ^ 
corresponds to flipping either a local maximum of ^ (a “corner” for the 
graph of ^) to a local minimum e vice versa. It is of course not positive 
recurrent, as the state space is inhnite and translation invariant for the 
dynamics, however, it is irreducible and recurrent. 

The reader can check that 0)^ is mapped onto Hw; by the transformation 
^ I— )■ where 


(6.9) vC(x):=C(x + l)-e(x) 

and that the image of the corner-flip dynamics (V^t)t>o is the simple exclu¬ 
sion process (see Figure 1). 

There is a natural order on the set n)y defined by 

(6.10) VxGZ2Ar, i{x)>^'{x), 


and we can construct a grand coupling for the Markov chain which preserves 
this order. 
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Eig. 1. The correspondence between the exclusion process and the corner-flip dynamics. 
A particle jump and its corner-flip counterpart are indicated by arrows. Note that this is 
not a one-to-one mapping as a particle configuration gives the height function only modulo 
translation. 

6.3. The graphical construction. We introduce in this section an order 
preserving grand-coupling on For (j £ (^^)t>o denotes the Markov 

chain with initial condition (j. We want to construct all the (^f)t>o on a 
same probability space in a way that 

(6.11) VC,C' e (C > C) => (Vt > > ef). 

Of course, there are several options for such a grand coupling. We want 
to choose one which is such that, eventually, the trajectories starting from 
different initial conditions coalesce almost surely (at a random time) 

(6.12) vc,c'G^^'iv,3rc,c'<oo>vt>rc,cL = 

Of course, we want the coalescing time to be as short as possible. To reach 
this aim, we make the different corner flips for different trajectories as inde¬ 
pendent as can be while still satisfying ( 6 . 11 ). 

Let us present the construction. The evolution of the is completely 

determined by auxiliary Poisson processes which we call clock processes. Set 

0 := {(x, z)\x £ Zjv and z £ 2 Z -|- (—1)^}. 

And set and to be two independent rate-one clock processes indexed 
by 0 (7J and 7j are two independent Poisson processes of intensity one of 
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each uj £ 0). The trajectory of given is given by the following 

construction: 

• is a cadlag, and does not jump until one of the clocks indexed by 
(x,^i(x)), X G I,2N rings. 

• If 7^^ ^ rings at time t and x is a local maximum for , then . 

rings at time t and x is a local minimum for , then = Cf- • 


6.4. Construction the initial condition for and Given x^Ga, 

we let the trajectory of the Markov chain with transitions rates (6.8) 
starting from initial condition 

X 

(6.13) ^0(a;):=^x(a;)- 

2=0 

Note that for all t > 0 we have 

(6.14) P[VC0g.]=P,^. 

Our idea is to construct another dynamic which starts with V^q dis¬ 
tributed like which coalesces with within time N'^^/logN. In fact, it 
turns out more practical to define not one but two dynamics and to 
couple with We let P denote the law of and we impose 

(6.15) p[VC,jG-]=P[V42G-]=i/“. 

Note that this implies for all t > 0 

(6.16) P[veiie.]=P[vei2E.]=i.«. 

We impose also the condition 

(6.17) < ^0 < eo^ 

and use the graphical coupling introduced in the previous section to con¬ 
struct the trajectory of (^t)t>o, * = 0,1,2. Hence, the order is conserved at 
all time 


(6.18) vt>o, 

Let us now explain our construction of the initial conditions. We start with 
Tfo distributed like and we will choose and such that 

(6.19) = %. 

We set for arbitrary rj G fljv, or ^ G 

y 




max 

x,y&N 


^ i]{z) - e~^^^sm{t) 

Z=X+1 


:=Ht,a{V0- 


( 6 . 20 ) 





We also set 
( 6 . 21 ) 
and 
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Uo ■= + VN\og\ogN)/2\ 


- no, 

, . Z = 1 

( 6 . 22 ) 

X 

+no. 

Z = 1 

The fact that (6.17) is satisfied follows from the definition of Qa and that of 
no- Note also that from Proposition 6.2 applied at t = 0, we have 

(6.23) F[no > 2VNloglogN] < (logiV)"^ 

To prove Proposition 4.2, it is sufficient to prove that and typically 
coalesce within a time N‘^^/\ogJ7. More precisely, we have the following. 


Proposition 6.5. For sufficiently large N, for all a < 2N for {f,l)t<o, 
{et)t>o, constructed as above, we have 


(6.24) 




1 

2 log log N' 


Proposition 6.5 is proved in Sections 6.5 and 6.6. 


Proof of Proposition 4.2. Let x Ga be fixed and consider the 
dynamics ff, i = 0,1,2 constructed as above. From (6.18), we have 

(6.25) => 

Recalling (6.14) and (6.16), we have for any t > 0 

(6.26) \\P^ - < IP[Ve? / ve'] < / (h- 

Hence, Proposition 6.5 implies the result. □ 

6.5. The randomly walking area. Let us set 

(6.27) A{t) = 

x&il 

The reader can check that A{t) is an integer. Because of (6.18), we remark 
that A{t) is always positive, and hence that and merge at time 

(6.28) r := inf{t > 0|H(t) = 0}. 
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As A(t) is an integer valued martingale which only makes ±1 jumps, it is 
to be a time changed symmetric nearest neighbor walk on In order to 
get a bound for 


F[T<t], 


we need to have a reasonable control over the time change, that is, the 
jump rate of A{t). It depends on the particular configuration the 

system sits on: it is given by the number of places where corners can flip 
independently for and More precisely, set 


(6.29) 


Ui{t) := {x € has a local extremum at x and 


3ye{x-i,x,x + 


The jump rate of A{t) is given by 


(6.30) u{t):=i^Ui{t) + #U2{t). 

For t < u{t)dt, let us define 


(6.31) 


J{t) := inf< s 


J n(n)dn>t|. 


By construction, the process {Xt)t>o defined by 

(6.32) Xt-.= A{J{t)) 


is a continuous time random walk on Z_|_ which jumps up and down with 
rate 1/2. From the definition, we have 

(6.33) Ao = A(0) := NUq. 


Note that from Proposition 6.2, and the definition of TLq we have 

(6.34) P[A(0) > 2N^/'^loglogN] < (logA)"^ 

To estimate r, we have to control the evolution of Xt (using standard 
properties of the random walk) and that of u{t) (using the properties of 
proved in Section 6.1). 


6.6. Multiscale analysis. To have the best possible control on u{t), we 
need to perform a multi-scale analysis. We construct a sequence of interme¬ 
diate stopping time (rj)i>o as follows: 

(6.35) Ti := inf{t > 0|A(t) < 

We set T_i := 0 for convenience. We are interested in Ti for i G {0,..., 
|'(log 2 A)/2]} where log 2 (-) := log(-)/log(2) denotes the logarithm in base 
2. To bound the value of r, we bound the value of each Ar^ = — rj_i for 

* < r (log 2 X)/2~\ and that of r — jv) / 2 ] • The way to do this is: 
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(i) First, we prove a bound for the analog of the Ar* for the process Xt 
defined in (6.32). 

(ii) Second, we prove a bound for u{t) which is valid in the interval 


For step (i), let us define 


(6.36) 


%:= 


Too- = 



It follows from standard properties of the random walk and from (6.34) that 
we have the following. 


Lemma 6.6. We have the following estimates: 

P[3i G {0,..., r(log2iV)/2l},7I > 3-'iV3(logiV)V4] < (logiV)-Vio, 

(6.37) 

P[Too > N^ilogN)^/^] < (logiV)-^/^o. 

For more details, we refer to the proof of [9], Lemma 6.1. 

Step (ii) is more delicate, because we cannot get a good bound on u which 
is uniform in time. For instance, we need to prove that most of the time u{t) 
is of order N but we know that just before r we have u{t) = 4. Hence, we will 
prove a different bound for each value of i. The bound is valid most of the 
time, and we will need to check that the small fraction of time during which 
it does not hold can be dealt with in the computations. Recalling (6.20), we 
set 

(6.38) m) := max(77„,i(ei^) + 

We notice that from the definition 

(6.39) - il{x)) < nit). 

Using this information, we can get the following control on u [recall (6.4)]: 
Lemma 6.7. If gT, we have 

(6.40) „(t)>in,in(]jV,4|i]). 

The proof is identical to the one of [9], Lemma 6.3. Note that thanks to 
Proposition (6.3) and our assumption a < 2iV“^/'^, the inequality (6.40) is 
valid up to time N'^ (which is much more than what we need) with high 
probability. To make this bound on u useful, we need to show that most of 
the time 'H{t) is not too large. 
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Lemma 6 .8. For any T > 0 , 


(6.41) P 


L4o 


^{H{t)>FNloglogN} 


-4 


< (logiV) 


-1 


Proof. It follows from (6.2) that for N sufficiently large, for any t > 0 


(6.42) F[n{t) > ViVlogloglV] < (logiV)-^ 

Then the result follows by using the Markov property for the integrated 
inequality. □ 


Proof of Proposition 6.5. Set 

(6.43) i3:=|^ l{w(i)>^/iviogiogA}‘^^ - 

C-.= {%< 3-'iV3(logAr)V4} n {Too < N\logN)^/^}, 

where 

(6.44) T := N'^y/logN. 

We assume also that N is large enough so that from Proposition 6.3 and 
Lemmas 6.6 and 6.8 we have 

(6.45) P[^n.BnC]>l-(21oglogiV)-^ 

Hence, the results follows if we can prove that 

(6.46) {AnBnC}c{T<T}. 

We split the proof of (6.46) in two statements. We want to show first that 
on the event Ar\B<^C 

(6.47) T - rpog2iv/2l < (logA^)^/'^iV^, 
and then that 

(6.48) VzG{0,...,r(log2iV)/2l}, {n - n-i) < {i + l)-^N\\ogNfl\ 
These inequalities combined give 

K 

(6.49) r < (logiV)i/^iV2 + ^(i + l)-^N^{logNy/^ < N^^logN. 

i=0 

Note that (6.47) is an immediate consequence of C as 

(6.50) Too= f u{t)dt>T — TK- 

Jtk 
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(6.51) 


Let us turn to (6.48). Let us assume that the statement is false and set 
io := minji e {0 ,..., 

r(log2iV)/2l}|(Ti -Ti_i) > (i + l)-2iV2(logiV)l/3}. 

The definition of io implies that 

(6.52) Ti,.i + (io + l)-2iV2(logiV)V3 < T. 

From B, we have [using (6.52) to obtain the second inequality] 

1 


'-{H{t)<VNloglogN} 


(6.53) 


rn^-l + iio+l) 2Af2(logAf)l/3 

- / '^{H{t)<VNlog\ogN}^^ 

■^Uo-l 

= (io + l)-2iV2(logAr)V3 

rT,o_i + (io+l)-=N2(logAf)i/3 


^{W(t)>\/AfloglogAf} 


-3 


> (io + ly^N^OogN)^/^ - N‘^{logN) 

>l(io + l)-2iV2(logiV)V3. 

For all t <TiQ, we have A(t) > , and thus using Lemma 6.7 and the 

assumption that A holds, 


u{t) > — min[ ■ 


Ait) 


(6.54) 


From (6.53), 


> 


8 \ ’ max('H(t), Ai^/2) _ 

j^3/22-io 

sVNlog log N ^{nt)<ANlog\ogN} ■ 


(6.55) 


r^o iV 2 “*o r*o 

^ - 8(loglogiV) l{«(t)<^/NloglogN} di 




This brings a contradiction to C (if N is large enough) and ends the proof 
of (6.46). □ 
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7. Proof of Proposition 4.3. To prove the result, we will try to control 
the derivative in t of the total variation distance that we have to bound. 

Note that ||r'“ — ^ ||tv is always differentiable on the right. This 

comes from the fact that for each 77 G Otv, both i'i{vi) and are 

differentiable. With a small abuse of notation, we use dt to denote the right 
derivative. Our method to prove Proposition 4.3 relies on getting a bound 
on the derivative valid for all a < 1. More precisely, we want to prove 

(7.1) dt\\vt - ^ 

Indeed, once integrated this gives 

(7.2) sup||i/f - ^^*||rpy < 03 ( 0 ^ + a^iV^/^), 

t>o 


which is equivalent to our result. 

Let us first perform a simple computation to show that it is sufficient to 
prove (7.1) in the case t = 0. Using the triangular inequality and the fact 
that the semi-group shrinks the total-variation distance, we have for any 
positive e, 


(7.3) 


Wt+e-^ IItV 


< ||,y« _ 


I TV 


+ 


< II,"^^‘ll _|_ 

— iTt ^ IItV ^ 


11 ,,oe“^iv* oe“'^Jv(*+®) II 

iTe ^ IItV 

ae~^N^ Q,g-7vjv(t+e) II 

^ IItV 


Hence, 


(7.4) 


,Q;e iV^ 


o ||_.a _.ae ''iV'ii ^ -• ,, 


Q,g—.^iV(^ + £) I 


lTVle=0' 


Note that the right-hand side is simply equal to 


9sIk" -V 




|rpv|S = 0 


for ol = ae Hence, to prove (7.1) it is sufficient to show that for all a < 1 

(7.5) dtM - 11 TV 14=0 < + C'2a^iV-3/2. 

We let gf denote the density of and g°‘ that of z/“. Recall that we have 

gaa{r]) 


(7.6) 

where 

(7.7) 

We have 

(7.8) 


9°{v)■= 


gj^{e^<v)y 


0 ( 77 ):= ^ r 7 (x)sin(a:). 

XSZ 2 JV 


9, 


'tm -4" 


.ae 


|TV = M7v|97(5f(47)-5"^ "(4?))|t=o|. 
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We compute the derivatives of gf and g'"® ^ {g) separately. We have 

(7.9) ^'^\g)\t=o = 0!>^Ng°'{g)[-a{g) + i^°‘{a{g))]. 

The other term requires more work, and we have 

(7.10) dtgnri)\t=o = Cg= ^ ff“(r?") - 5“(r?). 

x£Z2N 

Recall (6.9). We have 

(7.11) 5“(??"') - ff"(??) = ff“(??)[exp(-aVsin(x)V??(x)) - 1]. 
Performing a Taylor expansion of the exponential, we have 


(7.12) Cg:=g^ig) 


-ai ^ Vsm{x)Vg{x)j+—G{g,N) + R{g,N) 

XG^2N 


where {oP'/2)G{g,N) is the second term in the Taylor expansion 

(7.13) G{g,N):= ^ (Vihf(x))2(V77(x))^ 

xG'^2N 

and R{g, N) is the Taylor rest 

R(n, N) ;= Y, - 1 + aV^{x)Vg{x) 

X£.’L2N 

(7.14) ^ 

- Y(Vsin(i))^(V7((i))^). 

The first term in the RHS of (7.12) can be simplified using summation by 
part and the fact that sin is an eigenfunction of A. We have 

(7.15) Vsin(x)V7/(x) = — /S.sm.{x)g{x) = \Na{g)- 

X^'L2N X^’L2N 

According to Taylor’s formula, one has for all a < 1, for an adequate choice 
of constant Ci 

p2Q: 3 _ 

(7.16) |R(r?,iV)|<- V |Vsin(x)Vr/(x)|^ <Cia^iV-2, 

6 

xGZ2N 

where in the last inequality we have used that |Vr/(x)| < 2 and that 


(7.17) 


|Vsin(a;)| = 2 


sm 


vr \ f TTX vr 


TT 

< — . 

- N 
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Combining (7.8) with (7.9) and (7.15), we obtain 



(7.18) 



To estimate the first term, we note that as 


(7.19) 


/^iv(5t(5f(r/)-ff“'^""‘(r?))k=o) = 0. 


we have 



Hence, from (7.16) 


(7.21) 


Va Riv, N) - aAArz^"(a(r?)) + —i/„(G(r/, N)) 

= i^a\R{v,N) - u^{R{ 7],N))\< Cia^N-^. 


To estimate the second term, we use Proposition B.l. The reader can check 
that the Lipshitz norm of G [cf. (B.l)] of G satisfies 


(7.22) ||G(-,iV)||iip<87r2iV-2 


and hence that for an adequate choice of 6*2 > 0 


(7.23) v‘^\G{r],N) - i/"(G(?7,iV))| < 


This completes the proof of (7.5). 

APPENDIX A: PROOF OF TECHNICAL STATEMENTS ON i/". 

A.l. Proof of Proposition 6.1. Note that if was replaced by the uni¬ 
form measure on {—1, l}^2jv (without the constraint of having N particles) 
then would be a product of independent Bernoulli, and the statement 
would be trivial to prove. 

What we have to control is that the constraint on the number of particles 
does not affect the mean too much. To do so, we perform an expansion of 
the partition function according to the value of ri{x) to show that the ratio 
of the partition function restricted to the event rj{x) = -|-1 and r]{x) = — 1 , 
respectively, is close to exp(2asin(x)). To this purpose, we introduce the 








CUTOFF PROFILE FOR EXCLUSION ON THE CIRCLE 


27 


quantity 

(A.l) 

We have 


|ry(a;) = —1) 


i/“(r/(x)) 

= h'°‘{r]{x) = +1) — i'°‘{ri{x) = —1) 

(A-2) _ _ 

= + 1 ) — = — 1 ) 

~ //Ar(e“^-« 2 iv |,y(x) = +1) + |r?(x) = -1) 

_ e2“^(^)z(x) - 1 

“ g2«sin(x)2-(a;) + i' 

Hence, what we must check to prove (6.1) is that Z{x) is very close to one. 
Now note that we can obtain a coupling of fiNi-\ilix) = —1) and /i 7 v(-|?/(ic) = 
+ 1) in the following manner: take rj^ with distribution fij\f{-\r]{x) = —1), 
choose y uniformly at random (and independent of in {z\ri^{z) = +1} 
and let rf be obtained from by exchanging the value at x and y (which 
are +1 and —1, resp.). A consequence of this coupling is that 


Z{x) 

(A.3) 




and hence we can deduce from it 


Z(x) = f“(T ^ i±|Me-2«.ln(.)|,,(^) = _l) 


(A.4) 


y&iN 

= 1 + z/“ ( — V 

' N ^ 2 

y&2N 


1 ~L yiy^ ^^—2asin{y) _ 


1 ) rj{x) = —1 


Note that with this expression it is not hard to check that \Z{x) — 1 | < 
e^“ — 1. However, to get a sharper estimate, we must have a good control on 
z/“(T/(?/)|r/(x) = —1). We obtain it by pushing the expansion one step further. 
We set 


/i 7 v(e“^"'^* 2 Ar\U.!'}®‘°^^^|?^(x) = +1,77(2/) = —1) 

+7v(e'^^"’^*2iv\U.!/}®‘"^^^|?7(x) = —1,77(27) = —1) 
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Similar to (A.2), we obtain that 


(A.6) 


= - 1 ) = 




N e2“"“(2^)Z'(x,y) - 1 
iV _ 1 ^2c^{y)1i(x^ y) + 1 ' 

Like for (A.4), we have an alternative expression for Z' 

Z'{x,y) 

(A.7) 

Hence, we have 

(A.8) |Z'(a:,2/)-l|<e2“-l, 

and from (A.6), we deduce that for some positive constant Ci 

(A.9) ' .. ' ^ 

Hence, we have 


= -l,77(y) = -1 . 


\^°‘{v{y)\y{x) = -1)1 <Ci[— + a 


\Z{x)-l\ < 


2N 


(A.IO) 


(e-2asin(;/) _ 

y&2N\{x} 


1 




yGZ2N\{x} 


Performing a Taylor expansion up to the second order in ol we obtain (recall 

a < 1) 


(A.ll) 


(g-2Q;sm(y) _ 


< 2a|sin(x)| + 


eNor 


yeZ2jv\{a:} 

The second term in the RHS of (A.IO) can be bounded by 


(A.12) 

Hence, we obtain 

(A,13) 


C.l4 + a)(e“-l), 


\Z{x)-l\<C2{a^ + N-^). 


And then the result can easily be deduced from (A.2). 
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A.2. Proof of Proposition 6.2. The result follows from the combination 
of Proposition 5.1 which controls the fluctuation around the mean value 
vP°{x,t) given an initial condition 770 find the following statement, that the 
mean itself vP°{x,t) does not fluctuate too much if rjQ has distribution z/“. 

Lemma A.l. There exists a constant c such that for all N > 0, for all 
|q:| < and t>0 


(A. 14 ) i/" 

3 x, 7 / G Z 2 Ar, 

y 

y~~] {u^{x,t) — ae~^^^sm{z)) 

> sVn 



Z=X+1 



Proof. It is in fact sufficient to prove (A.14) for t = 0, because 

y y 

(A.15) max t) — ae“'^^*sin(z)) < max (r/o(- 2 ) — asin(z)). 

x^y ^ x^y 

z=x-hl z=x-hl 

Indeed, if one sets v{x,t) to be the solution of the discrete-heat equation on 
Z 2 N with initial condition 

X 

Vo{x) := — asin(2:), 

Z = 1 

then (A.15) can be reformulated as 


(A.16) max[7;(t, y) — v{t, x)] < max[uo(7/) — uo(x)] 


which is obviously true by contractivity of the heat equation. Note that at 
the cost of losing a factor in the constant c, we can restrict ourselves to 
proving that 


(A.17) 


By e Z2Ar, 




z=l 


>4sVn 


< 2e- 


We have used Proposition 6.1 and the assumption on a to replace asin(z) 
by v°^{r][z)). Let us introduce notation for the sum 

(A.18) Sa,^y-.= ^{r}{z)-v^{y{z))). 

Z=1 

We also set p := [log 2 -|- 1. For s > 0, we set 

J(s) := {3q e {1,... ,p},3y e {1,..., 


\S2i{y-l),2iy\ > (|)^ ‘^S'/N}. 


(A.19) 
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By a simple dichotomy argument (see the proof of Proposition 4.1 in [9]), 
we have 


(A.20) <3?/EZ2a/', 


y 

E 

z=l 


{r]{z) — 1/2 — asin( 2 :)) 


C J(s). 


For y and p fixed, S 2 <i(y-i)^ 2 iy{'n) is a function which depends on 2”? coordi¬ 
nates and whose Lipshitz norm is smaller than 2. Hence, by Proposition B.l, 
we have 


(A.21) 


< 2exp(-Ci(^)P-''2-‘'s2iV) 

< 2exp(-2C'i(|)^"'^s^). 


Hence, by a union bound, for an appropriate choice of constant C 2 and 
for all s > 0, we have 


(A.22) z/“CH(s)) < 2^2?’"''exp( -2Ci 

q=l 


p-q 


< 2exp(—C 2 S^). 


□ 


A.3. Proof of Proposition 6.3. Set 

(A.23) £ := {?7 e nAr|3(x,y) G Z|^,#[x,y] > ^ j{x,y,r]) > 

First, we notice that from the proof of [8], Lemma 6.2, there exists a constant 
Cl > 0 such that 

(A,24) 


Recall that zz" the law of rjt starting from distribution zz“. We have by the 
Cauchy-Schwarz inequality 


(A.25) 


{i.n£)f<m{£)m 



Note that the term /^Az[(g^)^] is decreasing in t, because the semi-group of 
the Markov chain contracts the I 2 norm. For t = 0, we have 


(A.26) 


^J-N 





Using Proposition B.l to have Gaussian concentration for a{9), we have for 
N sufficiently large: 


fJ-N 



< exp(lOOA^a^). 


(A.27) 
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Hence, we can conclude that there exists constant C 2 and Cs such that if 
a < C 2 N~^/^ for any t we have 

(A.28) uf{B) = [r]t e£]< e-^3iVi/4^ 

Now we have to move from this result to a result for all t < . Note that 

starting from rj ^ one needs at least transitions in order to jump 

out of £. Hence, using union bound 


(A.29) 


sEi 

i=0 


7V5 

[Vi/N^ c f ^ p 


2=1 


(^t)te[(i-l)/N2,i/Ar2] 


performs more than transitions 


The hrst term is smaller than cf. (A.28). As for the second one, it 

is not difficult to check that the rate at which transitions occur in the chain 
is bounded by 2N, and thus that for any i 

(A.30)P[(77f)jg[(j_i)/7v2,i/iv2] performs more than transitions] < e~^, 

provided N is large enough. 


APPENDIX B: CONCENTRATION EOR LIPSCHITZ EUNCTION OP 

PARTICLE SYSTEMS. 


Given / : {0, —)• M, one sets ||/||iip to be the Lipschitz norm of / for 

the Hamming distance 


(B.l) 


||/||lip:= max 


\f{'n)- f{i)\ 


Yjx&2N {^)} 


Proposition B.l. For any /{-l,l}^2iv we have 
(B.2) W(l/ - M«(/)l > ») < 2exp(- _'’p|| j:|| 2 ^ 


If the function f only depends on {rjx)xeA where A is fixed a subset of 'L 2 N 
of cardinal k we have 

(B.3) FN{\f - hN{f)\>s) 

The result remains valid if hn is replaced by a measure v whose density 
with respect to pN is of the form 




(B.4) 
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where g is an arbitrary function on 'Z 2 N ■ 

Proof. We can without loss of generality assume that ||/||iip = 1. Now, 
we introduce the martingale dehned by 

(B.5) Mi{g) := 

We are going to check that the increments of M are bounded, that is, 

(B. 6 ) Vie{0,...,2A^-2}, \Mi+i-Mi\<2 

and the proposition is then simply a consequence of Azuma’s concentration 
inequality [ 1 ]. 

To check (B. 6 ), we need to show that for any realization one 

can couple 77 ^ and 77 ^ with law 

7^1 ■=u{-\{g{x))l^^,ri{i + l) = l), 

(B.7) 

:=i 7 (-|( 77 (x))*^^i, 77(7 + 1 ) = -!) 

in a way that ( 77 ^ — rf‘){x) has only two discrepancies, one at z + 1 and 
another one in where 77 ^( 3 :) = 1 — 77 ^ (x) = 0 . 

Note that vi and V 2 can be considered as a measure on { — 1,1} ^ {i + 
2 ,..., 2 A^}, one which is concentrated on the set of configurations with k := 
N — ~ 1 particles, and the other on the set of configuration with 

k + 1 particles. What one can do is to first draw 77 ^ according to vi , and then 
add a 1 chosen at random to the configuration to obtain rf. One 77 ^ is given, 
and we choose at random a site X in {x £ {i + 2,..., 2 A^}| 77 ^(x) = — 1} with 
distribution 

p9{^) 

On can check that 77 ^ dehned by 

(B.9) 77^(x) := 77^(x) + l{x=x} - ^x=i+i}, 

has distribution 122 ■ 

For the case where / depends only on 771 ^, we can consider a fc-step 
martingale which unveils at each step the state r]{x) of one x £ A. □ 

Acknowledgments. The author is very much indebted to Milton Jara 
who suggested to him the question of cutoff prohle, and pointed out that 
the relaxation of the hrst Fourier coefficient of 77 should be similar to an 
Ornstein-Uhlenbeck process starting far from equilibrium. 
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